Photonic Boson Sampling in a Tunable Circuit 
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Quantum computers are unnecessary for exponentially-efficient computation or simulation if the 
Extended Church- Turing thesis — a foundational tenet of computer science — is correct. The thesis 
would be directly contradicted by a physical device that efficiently performs a task believed to be 
intractable for classical computers. Such a task is BosonSamplinG: obtaining a distribution of 
n bosons scattered by some linear-optical unitary process. Here we test the central premise of 
BosonSamplinG, experimentally verifying that the amplitudes of 3-photon scattering processes are 
given by the permanents of submatrices generated from a unitary describing a 6-mode integrated 
optical circuit. We find the protocol to be robust, working even with the unavoidable effects of 
photon loss, non-ideal sources, and imperfect detection. Strong evidence against the Extended- 
Church- Turing thesis will come from scaling to large numbers of photons, which is a much simpler 
task than building a universal quantum computer. 



Quantum computation has attracted much attention be- 
cause of the promise of new computational and scien- 
tific capabilities. The most famous quantum algorithm 
is Shor's factoring algorithm [T|, which if realised will 
efficiently factor large composite numbers into their con- 
stituent primes, a task whose presumed difficulty is at 
the basis of the majority of today's public-key encryp- 
tion schemes. What is not widely appreciated is that the 
very existence of Shor's algorithm poses a fundamental 
trilemma: respectively, at least one of the foundational 
tenets of physics, mathematics, or computer science, is 
untrue [2]. 

Shor's algorithm states that efficient factoring can be 
done on a quantum computer, which is thought to be a 
realistic physical device. It may be that a scalable quan- 
tum computer is not realistic, if for example quantum 
mechanics breaks down for large numbers of qubits [3_ . If, 
however a quantum computer is a realistic physical device 
at all scales, then the Extended Church- Turing thesis — 
that any computational function on a realistic physical 
device can be efficiently computed on a probabilistic Tur- 
ing Machine — means that a classical, efficient, factoring 
algorithm exists. Such an algorithm, long sought-after, 
would enable us to break public-key cryptosystems like 
RSA. A third possibility is that the Extended Church- 
Turing thesis itself is wrong. 

How do we answer this trilemma? As yet there is no ev- 
idence that quantum mechanics doesn't apply for large- 
scale quantum computers — that will need to be tested 
directly via experiment — and there is no efficient classi- 
cal factoring algorithm or mathematical proof of its im- 
possibility. This leaves examining the validity of the Ex- 
tended Church- Turing thesis. One approach would be to 
take a task that was believed, on strong evidence, to be 



intractable for classical computers, and to build a physi- 
cal device that performs the task efficiently. This would 
directly contradict the Extended Church- Turing thesis. 

One such task is BosonSampling: obtaining a rep- 
resentative sample-distribution for n bosons scattered by 
some linear-optical unitary process, given by a mxm ma- 
trix [7 [H This task likely becomes intractable for 
large n, for reasons related to the fact that the ampli- 
tudes of n-boson processes are given by the permanents 
of nxn sub-matrices of U, and calculating the permanent 
is a so-called '#P-complete' problem |5] — a complexity 
class above even 'NP-complete'. Note that BosonSam- 
plinG itself is not thought to be #P-complete: the abil- 
ity to solve it lets us sample random matrices most of 
which have large permanents, but probably does not let 
us estimate the permanent of a particular, given matrix. 
However, by using the fact that the permanent is t^P- 
complete, Ref. [4] showed that even for the 'easier' task 
of BosonSamplinG, any fast classical algorithm would 
lead to drastic consequences in classical computational 
complexity theory, notably collapse of the 'polynomial 
hierarchy'. (A similar task, instantaneous quantum poly- 
time, IQP, was recently proposed in Ref.[7|: experimen- 
tally it appears far more difficult to achieve). 

Here we test the central premise of BosonSamplinG 
in practice, experimentally verifying that the amplitudes 
of n=3 photon scattering events are given by the per- 
manents of nxn sub-matrices of the unitary operator U 
describing the photonic quantum computer. We find the 
protocol to be robust, working even with the unavoid- 
able real-world effects of photon loss, non-ideal photon 
sources, and imperfect photon detection. We make use 
of a novel method for experimental characterisation of 
near-unitary evolutions which is both fast and accurate 



2 




FIG. 1: Experimental scheme for BosonSampling. (a) 

In BosonSampling both Alice and Bob must find the out- 
put distribution from some unitary, U , for multi-boson in- 
puts; Alice & Bob respectively possess classical and quantum 
resources, (b) The equivalent circuit diagram of our unitary. 
Three spatial mode inputs each encode two orthogonal polar- 
isation modes which can be arbitrarily combined by the uni- 
taries ai, a2 and as. A 3x3 unitary evolution, u^'^', interferes 
all three modes of the same polarisation, these are recombined 
at the output by 6i, 62 and 63. Note the swap gates between 
modes 2 and 5, which reflect that only modes of the same po- 
larisation interact in u'^'. (c) Experimental schematic. Two 
pairs of single photons are produced via spontaneous paramet- 
ric downconversion in a nonlinear crystal (BBO), driven by a 
frequency doubled (SHG) femtosecond laser (Ti:S) see section 
I of Supplementary Material. Up to three of the four photons 
are injected into the BosonSampling circuit; photons 1 and 
3 can be delayed or advanced with respect to photon 2 by 
An, Ars respectively; the fourth photon serves as a trigger. 
The BosonSampling circuit is constructed from a combina- 
tion of 3x3 and 2x2 integrated optical fibre beam-splitters. 
The local unitaries, a\...bs, are implemented with polarisation 
controllers (POL); ^(3) is implemented by a biased 3x3 non- 
polarising fritter (FBS), its outputs are mapped to 6 spatial 
modes by three polarising fibre beam-splitters (PBS). The 
outputs of the final 6x6 circuit are coupled to single photon 
avalanche diodes (APDs) whose signals are processed by a 
counting logic based on a field-programmable gate-array cir- 
cuit. 



Imagine a race between two participants, Alice, who 
only possesses classical resources, and Bob, who pos- 
sesses quantum resources. They are given some uni- 
tary, U , and agree on a specific n-boson input config- 
uration. Alice calculates an output sample-distribution 
with a classical computer; Bob either builds — or pro- 
grams an existing — linear-photonic network, sending n 
single-photons through it and obtaining his sample by 
measuring the output distribution. The race ends when 
both return samples from the distribution. In principle 



the validity of Alice and Bob's samples can only be estab- 
lished by individually comparing them against knowledge 
of the full distribution; in practice they can check for con- 
sistency by seeing if their distributions compare to within 
error. The winner is whoever returns a sample fastest. 

As n becomes large, it is conjectured that Bob will al- 
ways win, since Alice's computation times increases expo- 
nentially, whereas Bob's experimental time does not. It 
becomes intractable to verify Bob's output against Al- 
ice's, and — unlike the case of Shor's algorithm — there 
is no known efficient algorithm to verify the result @]. 
Importantly, however, one can take a large instance — 
large enough for verification via a classical computer — 
and show that Bob's quantum computer solves the prob- 
lem much faster, thereby strongly suggesting that the 
same behaviour will continue for larger systems, casting 
serious doubt on the Extended Church- Turing Thesis. 

In a fair race. Bob will verify that his device actually 
implements the target unitary. An alternative fair ver- 
sion of the race is to give both Alice and Bob the same 
physical device — instead of a mathematical description — 
and have Alice characterise it before she predicts output 
samples via classical computation. Alice can use a char- 
acterisation method that neither requires nonclassical re- 
sources nor adds to the complexity of the task [5] . 

We performed our BosonSampling demonstration 
in an optical network with m—Q input and output 
modes, and n=2 and ?t.=3 photon inputs. We im- 
plement a randomly-chosen unitary which is fully- 
connected, i.e. every input is distributed to every 
output. The experimental setup is shown in Fig. [T] 
The 6-input X 6-output modes of the unitary U are rep- 
resented by two orthogonal polarisations in 3x3 spa- 
tial modes of a fused-fibre-beamsplitter, an intrinsi- 
cally stable and low-loss device. The mode mapping 
is {l,...,6M|i7)i,|y)i,|iI)2,|F)2,|i/)3,|^)3}, where 
\H)i is the horizontally polarised mode for spatial mode 
1. We use polarisation controllers at the inputs and out- 
puts to modify the unitary, see the equivalent circuit di- 
agram in Fig[Ij3). 

Previously, quantum photonic circuits have been char- 
acterised via quantum process tomography [9], requir- 
ing an exponentially-increasing number of measurements. 
A recent theoretical proposal — yet to be experimentally 
realised — suggests using a linear number of non-classical 
interferences |10j . however the output coincidence signal 
falls off quadratically with circuit size. Here we use an 
efficient, method for this challenging task which requires 
just (2m— 1) measurement combinations of single- and 
dual-mode coherent-state inputs to the mxm network. 
The photonic network is described by a mapping between 
the input, a|, and output, aj, creation operators, 

m 

a]=^t/,^,al (1) 
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FIG. 2: Two-photon BosonSampling. (a) Outline of the technique Alice uses to predict the visibility from the unitary 
evolution U. (We present the fully characterised unitary U in the section II of the Supplementary Material as well as an 
example permanent calculation). Alice's prediction is shown by the blue-line envelope at far right; the light-blue box represents 
the uncertainty, obtained by 10 separate characterisations of the unitary, (b) Sample two-photon quantum interferences: the 
five output combinations {l,m} for the input configuration of {1,5}. Errors are smaller than marker size and the solid blue 
lines are Gaussian fits used to calculate the visibility from Eq. [s] (c) Experimental data showing Alice's (predicted) and Bob's 
(measured) visibilities for two-photon quantum interference. The input configurations are shown in the top left of each panel; 
the output modes are labelled at the bottom of the plot. The solid blue-line envelopes are Alice's predictions for visibility 
based on her measurement of U; the orange bars are Bob's measured visibilities; the yellow circles are the visibility predictions 
if coherent input-states were used instead of two-photon inputs. Errors are given by the light-blue and dark-red boxes at the 
extrema of each data set. 



where Uij—ri 



' . In brief, Alice can determine the 
moduli Vij, by preparing single mode coherent states 
into modes l...m and measuring the output intensi- 
ties, which are proportional to the occupation numbers 

Tin . ..Tlij . . - Tiim , 

fij = y/rhj- (2) 

Next she determines the phases Oij by preparing pairwise 
interferometrically stable dual-mode coherent states be- 
tween outputs 1 and all other outputs 2...m. Further de- 
tails, including a full theoretical derivation, can be found 
in ref. [5]. 

Having obtained U, Alice calculates the probability of 
bosonic scattering events in the following way [H [TT] . 
Given the input and output configurations S={si, Sm) 
and T={ti, ...,tm) with boson occupation numbers Si 
and tj respectively, she produces an nxm matrix Ut 
by taking tj copies of the j"^ column of U. Then, 
she forms the nxn submatrix Ust by taking Si copies 
of the i"^ row of Ut- The probability for the scat- 
tering event T, for indistinguishable input photons S, 
is given by P^=\Pei{UsT)\'^ ■ Conversely, the clas- 
sical scattering probabilities — when the input photons 
are distinguishable — is given by P^— Per {Ust), where 

I 1 2 

UsTij= \ UsTij I ■ 

Bob on the other hand experimentally prepares the n- 
photon Fock state \ti, ...,tm)- After injecting the desired 
input to the circuit, he determines the probability of the 
scattering event T by projecting onto its corresponding 
state using single-photon detectors connected to a coin- 
cidence counting logic, see Fig. 



We prepare near-single-photon Fock states via spon- 
taneous parametric downconversion in a nonlinear uni- 
axial crystal. Fig. see section I in Supplementary 
Material for detail. Once the photons pass through the 
network, they are detected by single-photon avalanche 
diodes. The BosonSampling protocol measures the 
frequency of output events, i.e. raw coincident photon 
counts. 

These however are strongly affected by differences in 
efhciency between photon counters. To remove this effect 
we measure the non-classical interference visibility, 



Vt = 



pC _ pQ 



(3) 



where and Pj? are the quantum and classical proba- 
bilities for the output configuration T measured for com- 
pletely indistinguishable and distinguishable photons re- 
spectively. Distinguishable statistics are obtained by in- 
troducing a temporal delay, Ar, between the input pho- 
tons. When all photons are delayed by significantly more 
than their respective coherence lengths, L, true two- 
photon quantum interference cannot occur. Figure [2j\) 
outlines the technique Alice uses to predict the visibility 
from the unitary evolution U. 

For n=2, the output T is monitored continuously as 
a function of the temporal delay between the two in- 
put photons; as is typically done with the well-known 
Hong-Ou-Mandel effect [12 . For n—3, however, the 
low four-photon count rates mean that such a measure- 
ment will be degraded due to inevitable optical misalign- 
ment and drift encountered over necessarily long exper- 
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FIG. 3: Three-photon BosonSampling. Experimental data showing Alice's (predicted) and Bob's (measured) visibilities 
for three-photon quantum interference. The input configurations are shown in the top left of each panel; the output modes are 
labelled at the bottom of the plot. The solid blue-line envelopes are Alice's predictions for visibility based on her measurement 
of U; the orange bars are Bob's measured visibilities; the yellow circles are the visibility predictions if coherent input-states 
were used instead of three-photon inputs. Errors are given by the light-blue and dark-red boxes at the extrema of each data 
set. 



imental runtimes. Therefore, for n=3 the probabilities 
Pj^ are obtained from just two measurement settings, 
P^(l)={-Aro„0,AToo} and {2)={At^,0, -At^}, 
where {ti,T2,T3} are the temporal delays of photons 1, 
2 and 3 with respect to photon 2, and Atoo^L/c. P^ 
is calculated as the average of these two probabilities to 
account for optical misalignment. Accordingly, are 
obtained with a single measurement of the output fre- 
quencies for completely indistinguishable photons, given 
by the delays {0,0,0}. 

As a first test of our system we compare Alice's 
and Bob's methods by injecting n=2 photons into the 
BosonSampling circuit. Figure[2|3) shows a representa- 
tive set of non-classical two-photon interference patterns 
as a function of the temporal delay. At. A more com- 
plete picture is given in Fig. [2p), here we obtain sample 
distributions for three combinations of input configura- 
tions and measure at all C(6, 2) (6 choose 2) possible 
output configurations. We compare Alice's and Bob's 
measurements using the average Li-norm distance per 
output configuration, £1 = ^^^^ '^t \^t~^t\^ here it 
is £1=0.027, showing excellent agreement between them. 

Next we show that if Alice uses her available 
resources — notably, coherent-states — to perform an anal- 
ogous experiment to Bob's she will not obtain the same 
results (see section III in Supplementary Material). Her 
coherent-state predictions — given by the yellow circles in 
Fig. ^p) — are clearly different to Bob's quantum mea- 
surements, with £1 =0.548. This large disagreement in- 
dicates that the quantum distribution is sufficiently dif- 
ferent from the classical distribution — so that this uni- 
tary isn't some special choice that is easy to sample 
classically — and highlights that Bob is accurately sam- 
pling from a highly nonclassical distribution. 

Figure [3] shows the results for 7i=3 BosonSampling. 
Here we find £i=0.122: we attribute the larger aver- 
age distance chiefly to the increased ratio of higher-order 
photon emissions in the three-photon experiment com- 



a) b) 1 |{1,3,5} 



50:50 BS 




Lxi Lxi Lxi in CD 
Lxf Lxf Lxf in in 
cnT CO in 



FIG. 4: Three-photon BosonSampling with colliding 
outputs, (a) Number resolution was achieved with a 50:50 
fibre beam-splitter at the output of mode 5 and an additional 
single-photon detector. Note that an imperfect splitting ra- 
tio for this fiber beam-splitter impedes only the effective ef- 
ficiency of our number resolving scheme [131 114j . (b) For an 
input configuration of {1,3,5}, and measuring two-photons 
in output mode 5, the solid blue-line envelope shows Alice's 
predictions; the green bars are Bob's measured visibilities; 
and the yellow circles are the visibility predictions if coherent 
input-states were used instead of three-photon inputs. Errors 
are given by the light-blue and green boxes at the extrema of 
each data set. 



pared with the two-photon case (see section IV in Supple- 
mentary Material). Again, Alice's efficiently-computed 
coherent-state predictions are clearly different to Bob's 
measurements, with £i=0.358. 

Having tested all possible 'non-colliding' output 
configurations — that is, 1-photon per output-mode — we 
also tested 'colliding' configurations with up to 2-photons 
per output-mode. This requires photon-number reso- 
lution [131 El J using the method shown in Fig. [4]A). 
Figure |4|3) shows once again good agreement between 
Alice's predicted and Bob's measured sample distribu- 
tions, £1=0.153, and a much larger distance between 
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Bob's measurements and Alice's coherent-state predic- 
tions, £1=0.995. 

Alice's Fock-state predictions and Bob's measurement 
results do not quite overlap to within error for n=2, and 
more notably disagree for n=3. This indicates the pres- 
ence of respectively small and moderate systematic dif- 
ferences between Alice's and Bob's methods for obtain- 
ing non-classical visibilities. This is as expected, since 
Alice's calculations are for indistinguishable Fock-state 
inputs, and Bob does not actually have these. It is well- 
known that the conditioned outputs from downconver- 
sion have higher-order terms, i.e. a small probability 
of producing more than one photon per mode (see sec- 
tion IV of the Supplementary Material), and are also 
spectrally-entangled, leading to a small degree of distin- 
guishability. The higher-order terms increase with source 
brightness, as Fig.[5]shows, this increases the distance be- 
tween Alice's Fock-state predictions and Bob's measure- 
ments, pushing the latter into a more classical regime. 
Thus when using downconversion, source brightness must 
be kept low, but since downconversion is probabilistic 
its brightness decreases exponentially with n: the best 
demonstrations to date have been n—8 at a rate of ~10~'^ 
Hz [T5]. Some gains can be made by increased spec- 
tral filtering, which decreases photon distinguishability. 
As the top and bottom data points in Fig. [5] show, this 
significantly moves Bob's measurements towards Alice's 
Fock-state predictions and away from the classical case. 

Strong evidence against the Church- Turing thesis 
will come from demonstrating BosonSampling with a 
larger-sized system where Bob's experimental sampling 
is far faster than Alice's calculation and where classical 
verification is still barely possible — according to [4], this 
regime is on the order of n=20 to n=30 photons in a 
network with with m^n modes. 

This is beyond current technologies, but rapid im- 
provements in efficient detection, low-loss circuitry and 
improved photon sources are highly promising. Highly- 
efficient, photon- number-resolving, detectors are now 
well-developed: e.g. transition edge sensors [TH] have 
high intrinsic efficiencies, >95%, recently enabling whole- 
system detection efficiencies with entangled-photon pairs 
of >60% fT7|. Integrated circuitry is also being rapidly 
developed. Nonclassical interference of photon pairs has 
been observed in periodically-bounded [TB] and linear [TS] 
waveguide arrays, up to m=21. Whilst not strictly nec- 
essary, reconfigurable integrated circuits [2Ql [21] could be 
used to sample from an ensemble of different unitaries. 
However, if they are to be used for high-n demonstra- 
tions, these circuits will have to be significantly improved 
to reduce their very high losses, e.g. incurred by coupling 
photon sources to waveguide arrays. Finally, triggered or 
heralded photon-sources will be required: Fig. [5] makes 
clear that downconversion alone is not a suitable photon 
source. Developing better photon sources is a focus of 
research worldwide, again with promising recent devel- 
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FIG. 5: Imperfect Fock states in BosonSampling. For 

an input configuration of {1, 3, 5}, increasing the downconver- 
sion pump power introduces larger proportions of higher-order 
photon numbers, red dashed-line and data: The Li-norm dis- 
tance between Alice's predictions for three-photon inputs and 
Bob's measured visibilities, black dotted-line and data: The 
I/i-norm distance between Alice's predictions for coherent- 
state inputs and Bob's measured visibilities. The lines are 
best fits to guide the eye. Data shown in black and red were 
taken using a spectral filter on the downconverted modes with 
FWHM of 5 nm. Decreasing the spectral filter to 2 nm sig- 
nificantly reduces the Li-norm distance between Alice's pre- 
dictions for Fock-state inputs and Bob's measured visibilities 
(green), while increasing the distance between Bob's measure- 
ments and those predicted using coherent state inputs (blue). 



opments, e.g. quantum-dot photon sources with produc- 
tion efficiencies of >40% 22 , orders-of-magnitude higher 
than downconversion. BosonSampling could also be 
demonstrated in other engineered quantum technologies, 
such as superconductors, neutral atoms, or, in an attrac- 
tive recent proposal, trapped ions: the phononic modes of 
the ions are the bosons; the ions act as the programmable 
beamsplitters of the phononic network 23J . 

An important open question remains as to the practical 
robustness of large implementations. Unlike the case of 
universal quantum computation, there are no known er- 
ror correction protocols for BosonSampling, or indeed 
any of the models of intermediate quantum computation, 
such as deterministic quantum computing with one qubit 
(DQCl) [531[5S], temporally unstructured quantum com- 
putation (IQP) [7], or permutational quantum comput- 
ing (PQC) [26]. These intermediate models have gar- 
nered much attention in recent years due both to the 
inherent questions they raise about quantum advantage 
in computing, and because some of them can efficiently 
solve problems believed to be classically intractable, e.g. 
DQCl has been applied in fields that range from knot 
theory [27] to quantum metrology [28] . Our experimental 
results are highly promising with regard to the robustness 
of BosonSampling, finding good agreement even with 
clearly imperfect experimental resources; equally heart- 
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ening, a recent theoretical study provides evidence that 
photonic BosonSampling retains its computational ad- 
vantage even in the presence of loss [29] . 
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Experimental Boson Sampling: 
Supplementary Material 



I. Four-photon source and BosonSampling circuit 

A mode-locked Ti:Sapphire {Coherent 900 HP) laser 
with a repetition rate of 76 MHz, 100 fs pulses and an 
average output power of ~3.8 W at 820 nni is frequency 
doubled in a 2 mm long bismuth borate (BiBO) nonlinear 
crystal to give ~1.5 W centred at 410 nm. Photon-pairs 
are produced via spontaneous parametric downconver- 
sion: the two pairs are produced by forward and back- 
ward passes of a 2 mm long beta-barium-borate (BBO) 
nonlinear crystal cut for type-I phase matching. Single 
photons pass through spectral filterers (FWHM 5 nm) 
before being coupled into single mode optical fibres. At 
100% pump power the forward and backward passes of 
the source produce approximately 290 kHz and 180 kHz 
two-fold coincidences, with a maximum four-fold rate of 
1.20 kHz. The discrepancy in two-fold coincidences be- 
tween forward and backward passes is due to differing 
focussing conditions in each case. 

Single photons are directed into the BosonSampling 
circuit using a combination of calcite beam-displacers 
and waveplates The forward pass of the downconversion 
source is used for the two-photon measurements and the 
source is run at 20% of the maximum pump power to re- 
duce the ratio of higher-order photon emission. With an 



average input coupling efficiency of 64% into all modes of 
the circuit, we obtain an average two-photon coincidence 
rate of 260 Hz across all modes. When injecting three 
photons into the circuit, the remaining fourth photon is 
sent directly to a single photon detector to serve as a trig- 
ger for its twin. Again, running at 20% maximum pump 
power we obtain an average four-fold coincidence rate of 
185 mHz. For all measurements in the main text we use 
an additional 2 nm filter at the input of the BosonSam- 
pling circuit, this reduces the effect of group- velocity 
mismatch between interfering photons from independent 
downconversion events. 

II. Measured unitary matrix and calculation of 
permanents 

In this section we present a sample measured unitary 
evolution and demonstrate how we generate sub-matrices 
whose permanents are used to calculate the probability 
amplitudes of bosonic scattering events. 

Changes in laboratory conditions cause slight changes 
to the in the optical fibres network we use as our Boson- 
Sampling circuit. For this reason we measure the evo- 
lution J/oxpj of the linear optical network after each of 
Bob's experimental runs, ensuring that Alice generates 
the most representative scattering probabilities. Below 
we give the unitaries that were measured after we com- 
pleting the two- and three-photon non-classical interfer- 
ence measurements respectively. 



7-r2photon 
^exp 



0.; 



'0.297 

,330 
0.388 
0.311 
0.396 
Vo.279 



0.325 - 
-0.302 

0.182 - 
-0.220 
-0.222 

0.322 - 



- O.OOOi 
' colli 

- 0.248i 

- 0.315i 

- 0.169i 

- 0.244i 



0.126- 
0.001 - 
-0.220- 
-0.169 - 
0.387- 



- O.OOOi 

- 0.503j 
-0.133i 

0.246i 
0.130j 



-0.101 - 0.239i 



0.500 + O.OOOi 
0.028 - 0.390j 
-0.212-|-0.204j 
0.190-l-0.157j 
-0.265-l-0.004j 
-0.051 - 0.400j 



0.430 + O.OOOi 
0.22H-0.118i 
-0.127- 0.386j 
-0.073 - 0.089j 
-0.103 -f0.202i 
-0.184 -f0.320i 



0.253 + O.OOOA 
-0.385 - 0.213i ' 

0.108 - O.OSli 
-0.227 4- 0.355i 

0.353 - 0.112i 
-0.217 + 0.074i/ 



7- rSphoton 



/0.334 
0.273 
0.420 
0.284 
0.340 

\0.324 



0.277 
-0.329 

0.140- 
-0.197 
-0.329 

0.344- 



O.OOOi 
0.051i 
0.242i 
0.367i 
0.049i 
■ 0.036i 



0.125- 
0.055- 
-0.191 - 
-0.194- 
0.328 ■ 
-0.114 



- O.OOOi 

- 0.478j 
-0.198j 

- 0.224i 

- 0.312i 

- O.lOli 



0.479 + O.OOOi 
0.021 - 0.121i 
-0.195 -f0.204i 
0.189 -(-0.190j 
-0.144 0.042 j 
-0.206 - 0.398i 



0.415 -(- O.OOOi 
0.197 -f0.128i 
-0.139 - 0.393i 
-0.072 - 0.106j 
-0.131 + 0.187* 
-0.111 + 0.351i 



0.237 + O.OOOA 
-0.345 - 0.253i 

0.113 - 0.085i 
-0.278 + 0.333i 

0.283 - 0.216i 
-0.098 + 0.208i/ 



For our example calculation we will input photons with 
the configuration S'=:(1,0, 1,0,0,0) and look at the out- 
put modes given by T=(0, 1, 0, 0, 1, 0). First we generate 
the sub-matrix Ut by selecting tj copies of the j"^ column 
of C/^i'p*'"*™ H E] to give. 



I 0.325 + 0.000i 
-0.302 - O.Olli 
0.182 -h0.248j 
-0.220 - 0.315j 
-0.222 - 0.169i 

V 0.322 + 0.244j 



0.430 -I- O.OOOiX 
0.221 -fO.llSi 
-0.127 - 0.386i 
-0.073 - 0.089i 
-0.103 + 0.202i 
-0.184 + 0.320i/ 
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Then we generate Ust by selecting st copies of the i^^ 
row of Ut to give 



0.325 
0.182 



O.OOOi 
0.248i 



0.430 + O.OOOi 
-0.127 - 0.386i 



(5) 



The probabihties of observing the even T for completely 
indistinguishable and distinguishable photons are given 
by and Pj?" respectively, 





|Per([/5T)|' 


(6) 




0.0017 




pC _ 


Per (Ust) 


(7) 




0.0349 




|2 

^ij 1 ■ 


The above probabilities 


give a 



where UsTij = \ Us' 

predicted non-classical interference visibility of V={P^ — 



P^)/pC^0.951, 



which compares well to Bob's observed 
value of 0.939 for this input/output configuration. 

The estimated errors for predicted visibilities given in 
the main text were calculated as the standard deviation 
of the predicted visibilities from 10 separate unitary char- 
acterisations. 

III. Calculation of visibility using coherent-state 
inputs 

We calculate interference visibilities by injecting n 
equal-amplitude coherent states with normalised elec- 
tric field amplitudes Ei—e'^^', into the input modes 
{i, m} of the linear optical network. The input vector 
E={Ei, Ei..., En) is transformed under U such that 
the electric field at output mode j is given by 



E, 



E,, 



(8) 



When the input coherent states overlap with a zero time 
delay the 2n-order correlation in electric field between 
detectors at output modes {j, m} is given by the phase 
averaged cross-correlation function [2TJ |30] , 



p(0) ^ 



1 

W) 



2-K 



2n rn 



■ dOi ...de„ 



(9) 



Conversely, if input coherent states are delayed by signif- 
icantly more than their coherence lengths no interference 
can occur between them, leaving the cross-correlation 
function as an incoherent sum of input fields at the out- 
put modes j, 



p(oo) 



(2^)" Jo 




de, ...de„ 



where E)j is the electric field at output mode j given 
the input at mode i. The interference visibility is then 
defined analogously to Eq.(3) in the main text, 



V = 



p{oo) _ p(0) 



(11) 



IV. Effects of non-ideal photon sources from 
dovirnconversion 

Spontaneous parametric downconversion, as the name 
implies, is a probabilistic phenomenon which outputs 
photon states of the form 



\tlj) cx |00) +r]\ll)+T]^ |22) -f- Tj^ |33) 



(12) 



where: \nani,) gives the photonic occupation numbers 
Ua and Uf) in the spatial modes a and b respectively |31| : 
the probability amplitude of creating one photon-pair per 
pump pulse is 77, which incorporates the nonlinear inter- 
action strength, interaction time of the pump field with 
the nonlinear crystal and the optical coupling efficiency 
of the down-converted modes, and r]<^l. Triggering us- 
ing one photon of a pair removes the vacuum component, 
but since 77* >0 there is a possibility that two or more pho- 
tons remain in the signal modes — these are the so-called 
higher-order terms. 

Higher-order terms can never be entirely removed from 
downconversion light with linear optics alone. Their 
detrimental effect can be reduced though by keeping 77 
at a necessary minimum, while keeping the single-photon 
rate constant through temporal 31 or spatial [32: mul- 
tiplexing of multiple downconversion sources. 



